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Abstract 



A negatively charged quasi-two dimensional exciton {X ) is solved ex- 

m ' 

■ actly numerically in the presence of a uniform perpendicular B-field. Various 

N£> 

quasi-two dimensional geometries are studied. The charge distribution of the 
X~ parallel to the B-field is found to be crucial in determining the stability 
ON ■ of the optically-active X and hence its photoluminescence (PL) signature. 



The theory provides a quantitative explanation of recent experimental results 



a 

*^ • obtained for a GaAs quantum well. Effects are found which cannot be de- 

scribed within a lowest Landau level approximation. 



PACS: 78.20.Ls 78.66.-w 73.20.Dx 



Excitonic effects are used to explain the majority of discrete lines appearing in the 
interband optical spectra obtained from semiconductors. However new lines have been 
seen recently which cannot be explained by excitons, either intrinsic or impurity-bound 
(see eg Ref. [|1| [0). It has been proposed that these discrete lines are associated with 
the formation of a charged electron-hole complex. The stability of charged excitons was 
first predicted for three-dimensional (3D) systems by Lampert |§. This was subsequently 
studied theoretically at zero B-field in both 3D || and two-dimensions (2D) [§. In both 
cases the charged exciton was found to be stable. The 2D charged exciton binding energies 
are about ten times larger than than those in 3D due to the extra confinement H thereby 
facilitating their experimental observation. Negatively charged excitons, X~ = 2e + lh, were 
first observed in CdTe/CdZnTe quantum wells || |7j and were later seen in GaAs/GaAlAs 
quantum wells | | | | [^0| . Positively charged excitons, X + = le + 2h, have 
also been seen [jl| . The charged excitons have been observed experimentally using polarized 
spectroscopy, in particular photoluminescence (PL) with a uniform magnetic field B applied 
perpendicular to the quantum well plane. X~ stability requires a two-dimensional electron 
gas (2DEG) of intermediate electron density which is typically near the metal-insulator 
transition, the reduced screening effect of the other electrons helping to stabilise the X~ [Q. 
In the low density limit the exciton (X) dominates the PL spectrum as expected. At high 
densities and low B-fields, a Fermi-edge-singularity (FES) is seen. The FES moves smoothly 
into the X~ line as the electron density is decreased ||] and develops into sharp X and 
X~ peaks with increasing B-field 0. These results suggest that a study of the X~ is crucial 
for an understanding of the optical properties of 2DEGs at finite B-fields, as was claimed 
recently in Ref. [|TJ. 

In this paper an accurate numerical solution is presented for a quasi-2D X~ consisting 
of two electrons and a hole situated in a weak in-plane confinement potential and a perpen- 
dicular uniform B-field. The X~ properties are studied over a wide range of experimentally 
applicable B-fields. Various relevant quasi-2D geometries are considered, each having a dif- 
ferent charge distribution parallel to the B-field, ie. perpendicular to the 2DEG. These 
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geometries are as follows: 

(i) electrons and hole are strictly 2D but move on separate parallel planes (biplanar geome- 
try), 

(ii) electrons and hole occupy the same plane but have their charge smeared in a rod-like 
distribution along the perpendicular direction (rod geometry), 

(hi) electrons and hole are strictly 2D and occupy the same plane (coplanar geometry). 
The stability of the optically-active X~ and hence its PL spectrum are found to vary dras- 
tically depending on the geometry used. It is shown that recent experimental PL spectra 
obtained from quantum wells can only be explained using the rod geometry and higher Lan- 
dau level mixing. A previous theoretical study only considered a coplanar geometry at a 



fixed B-field [11 



The Hamiltonian for the X problem is: 
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where m e = 0.067m and = 0.475m are the masses of the conduction electrons and 
heavy holes respectively for GaAs. In order to facilitate comparison with a quasi-2D electron 
gas the parabolic confinement potentials, huj e = 1.75meV and fiu^ = 0.25meV, are chosen 
to be weak compared with the cyclotron energies (h(uj c ) e = 17.2meV and h(u c )h = 2.4meV 
at B=10T). Also uj e /u>h = rrih/m e ensuring that the electron and hole wavefunctions have 
equal characteristic lengths, l 2 eh = fi/ {m e ^oJ e ,h) = I 2 where u 2 h = {{oj c ) 2 eh + 4a;g h ), as is the 
case in the limit of zero confinement. 

The interaction potentials are Coulombic. In order to obtain a tractable model it is 
assumed that the single particle wavefunctions separate: \l/ e = (fi e (z e )ip e (r e ) and = 
^hi^^^h ( r h)- A specific form for is chosen, thereby freezing out the z-motion and 
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yielding a quasi-two dimensional model. For the biplanar geometry: 

|0 e (z e )| 2 = 5(z e -O) ; \Mzh)\ 2 = 5(z h -d) 



(2) 



while for the rod geometry: 



- if — < z u < - 



lM**)r = < ' ( 3 ) 

otherwise 



The coplanar geometry is the limit of the previous two where d and L tend to zero. 

The coplanar and rod models retain a singularity in V e ^ at r e = rh which the biplanar model 

loses. The differences in the forms of |</>(z)| 2 lead to significantly different PL spectra. 

The effects of the B-field and the Coulomb interaction compete: the expansion parame- 
ters are given by a e ^jl where a £y h = th 2 /m e ,h& 2 is the Bohr radius for the electron and hole 
respectively. Taking B ~ 8T and e r = 12.53 yields (GaAs) a^jl = 0.2 and a e /l = 1.1, thereby 
rendering any perturbation approach unreliable. For this reason the X~ Hamiltonian is di- 
agonalized exactly numerically. Rotational symmetry about the B-field leads to conservation 
of total angular momentum, M. The symmetric gauge is hence chosen, A(r) = (B x r)/2, 
giving single particle states characterised by the Landau level (n) and azimuthal angular 
momentum (m). Angular momentum conservation makes the Hamiltonian block diagonal 
in this basis, each block being characterised by M. The spin and spatial components of the 
wavefunction factorise, the spin states of the electrons being either singlet or triplet. In 
order to preserve the correct electron antisymmetry the triplet (singlet) states have a spatial 
component which contains antisymmetric (symmetric) combinations of the single particle 
electron wavef unctions. 

The Hamiltonian is diagonalized in the basis of definite M, S and S z . Since (a e> h/l) > 1 
higher Landau levels must be included in the basis. The actual basis used has 6 hole Landau 
levels, 2 electron Landau levels and 20 angular momentum states per Landau level (cf. [|TT|] ). 
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A larger number of hole Landau levels are used because of the small heavy-hole cyclotron 
energy. A finite number of angular momentum states are used due to the weak confinement 
which removes the degeneracy of the Landau level. From comparison with experiment 
(discussed later in this paper) the authors believe that this basis captures the essential 
physics of the X~. The matrix elements are evaluated following a method outlined in W2 



In the coplanar limit (L, d — > 0) analytic results are obtained for the matrix elements [13 
llj. Away from this limit, numerical integration is employed. Two completely different 
methods for evaluating the matrix elements have been used to check for numerical error in 
the integration routines. For comparison purposes the exciton lines are also calculated in a 
manner identical to that used for the X~ . 

Both excitons (X) and X~ undergo recombination giving rise to a PL spectrum. It is 
assumed here, as in 0, that recombination is with the ±3/2 heavy holes. The oscillator 
strengths for the X and X~ PL lines are (0|L|X) and (e\L\X~) respectively, where L = 
J ■ip e (r)'iph(r)d 2 Y. This means that only excitons with orbital angular momentum M=0 will 
recombine. For both X and X~ the necessary change in the system's angular momentum 
is taken up by the atomic part of the wavefunction giving rise to two lines, a + and o~~, 
which are Zeeman split. In the limit that the basis is restricted to the lowest Landau level 
and setting V e h = V ee , there is a well-documented hidden symmetry |14| Hl5fl . This forces 



the triplet X~ (Xf) PL energy to be identical with that of the lowest exciton (X) line. In 
this paper the symmetry is broken by the weak confinement and more importantly by the 
inclusion of higher Landau levels. For the case of a biplanar z-geometry the interactions V ee 
and V e h are no longer equal, thus further breaking the hidden symmetry. The reliability of 
the calculation has been verified by reproducing known results as follows: (a) in the limit 
of zero confinement and using the rod geometry (V ee = V e h) the hidden symmetry result is 
recovered, (b) the B=0 confined 2-electron results of |Uj are reproduced, (c) results found 



in a previous X~ study [|TT| are obtained. 

The theoretical PL spectra resulting from the coplanar, biplanar and rod geometries are 
now discussed. When a strictly coplanar geometry is used both the optically-active singlet 



(X~) and triplet (Xf) states are stable with respect to the decompositon X~ — > X + e. 
This stability persists over a calculated range in B-field of at least 6 — > 14T. However, the 
agreement with recent experimental quantum well data (see Fig la) is not good since the 
calculated X — X~ splitting is too large (3.8meV at B=8T) and there is no strong increase 
of the PL energies with B. Using the biplanar geometry the PL spectrum is found to be 
very sensitive to changes in the plane separation, d. For an X~ with zero total angular 
momentum, M=0, the optically-active singlet and triplet species are stable when d = 0. 
Increasing d from zero causes the singlet (X~) to unbind (d ~ 0.5/) with the triplet (Xf) 
unbinding soon after (d ~ I). Similar unbinding occurs for the optically-active M^O X~ 
studied. Thus the biplanar geometry PL spectra differ strongly from the experimental 
quantum well PL spectra (Fig la). These results demonstrate that quasi-2D systems which 
approximate to a biplanar system (eg heterojunctions) with d > I are most unlikely to 
exhibit PL effects due to charged excitons. 

The experimental PL spectrum of Shields et al. for a GaAs quantum well || shown 
in Fig la is typical of those found in the literature. The exciton (X), singlet X~ (X~), 
and triplet X~ (Xf) are marked in accordance with the experimental assignments of Ref. 
H (for clarity, only the a~ lines are shown in all figures in this paper). The theoretical 
lines shown in Fig lb are calculated using the rod geometry with a rod length L =225A 
and total angular momentum M=0. For the range of B-fields shown in Fig 1, L is less 
than the width of the X~ in the plane thus validating the assumption that the X~ is 
quasi-2D. The L value is also consistent with the width of the envelope function in the 
experimental quantum well used by Shields et al 0. Figure 1 shows that the theoretical 
lines exhibit many of the quantitative features of the experimental spectrum. The energy 
of the calculated PL lines increases with B-field in a manner closely matching that of the 
experimental lines. This is a non-trivial agreement since the slope is not the result of simply 
adding the single-particle cyclotron energies. The X — X~ splitting (~ 1.4meV at 8T) is 
close to the experimental value (~ 1.7meV at 8T). With increasing B-field the triplet line 
(Xf) becomes more separated from the exciton line (X) hence Xf becomes increasingly 
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stable. However, Xj~ only becomes stable at a finite B-field (5T theoretically compared 
with 2.5T experimentally). The ordering of the PL lines in energy is very stable to changes 
in L. The quantitative theory discussed here thus allows a definite assignment of the PL 
lines which is in agreement with that proposed for the original experimental spectra [0. 
The slight differences between the theory and experiment in Fig 1 are to be expected since 
the model is necessarily a simple, quasi-2D representation of a complex 3D system. 

Recombination from charged excitons in the rod geometry with total angular momentum 
M<0 should also be considered as they give PL lines of similar energies to M=0 X~ . The 
PL spectra for a M=— 2 X~ at B=10T are calculated both with and without an in-plane 
confinement and are shown in Figs 2a and 2b respectively. As the in-plane confinement is 
removed the lowest triplet state has its oscillator strength drastically lowered to the benefit 
of the first excited triplet state (cf Figs 2a and 2b). The PL spectra in the zero confinement 
limit agree well with the experimental PL spectra in Ref. ||. The samples used in Ref. 
are reported to have a thick undoped spacer layer between the donor atoms and the well. 
This minimizes the spatial confinement in the plane of the well, giving mobilities which 
demonstrate a high degree of homogeneity and allowing comparison with calculations in the 
zero confinement limit. Fig 2c shows the PL spectrum of the M=0 X~ for B=10T with 
the weak in-plane confinement. As the in-plane confinement is removed (Fig 2d) the PL 
spectrum is almost unaffected. If the basis is restricted to the lowest Landau level, the 
hidden symmetry result (triplet and exciton PL lines coincident) is obtained when the in- 
plane confinement is removed (see inset in Fig 2). The fact that the splitting is retained 
in Figs 2b and d implies that higher Landau levels are crucially important in allowing the 
optically-active Xf wavefunctions to form configurations stable against the decomposition 
Xf — > X + e~: this breaks the hidden symmetry and causes the Xf PL line to be seen. 
In fact the energy levels of the X~ and exciton all change significantly when the basis is 
restricted to the lowest Landau level. 

To summarize, a quasi-two dimensional X~ has been solved exactly numerically in the 
presence of a uniform perpendicular magnetic field and a weak in-plane parabolic confine- 
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ment. The diagonalization basis contained higher Landau levels giving rise to results which 
differ significantly from those obtained within a lowest Landau level approximation. The 
quasi-2D coplanar, biplanar and rod geometries studied were found to yield significantly dif- 
ferent PL spectra. Of all the geometries considered, the rod model approximates closest to a 
quantum well and was able to produce good quantitative agreement with recent experimen- 
tal quantum well spectra. The biplanar model predicts an instability of optically-active X~ 
with increasing d: this is consistent with the apparent failure to observe charged excitons in 
experimental quasi-2D systems other than quantum wells. 

The authors would like to thank Darren Leonard for illuminating discussions and J. J. 
Palacios for a helpful comment. Funding was provided by EPSRC through a studentship 
(J.R.C.) and a Research Grant No. GR/K 15619 (N.F.J, and V.N.N.) . 
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Figure Captions 

Figure 1: Variation of X~ and exciton X PL spectra (<x~) with B-field. (a) Experimental; 
adapted from || for a 300 A GaAs quantum well, (b) Theory; zero of energy shifted 
to include the band gap. Ranges are chosen to highlight the agreement of trends in the 
experimental and calculated spectra. Parameters given in text. Lines are guides to the eye. 



Figure 2: PL spectrum (cr~) at B=10T, L =225A for (a) M=— 2, weak in-plane confinement 
(b) M=— 2, zero in-plane confinement limit (c) M=0, weak in-plane confinement (d) M=0, 
zero in-plane confinement limit. Parameters given in text. Inset shows Xf and X lines 
calculated in the lowest Landau level (LLL) with zero in-plane confinement. A gaussian 
line-broadening has been introduced throughout. 



10 



1525.0 



QJ 

B 

>» 

u 

QJ 

a 

w 

QJ 
S3 
w 

S3 

a 

s 

o 



1524.0 



1523.0 



1522.0 



1521.0 



1520.0 



• Exciton X 

♦ Triplet X" 
■ Singlet X 



• ♦ 
• ♦ 

• ♦ 
• ♦ 



• ♦ 



0.0 



2.0 



4.0 



1522.0 



1521.0 



1520.0 



2.0 



(a) 



6.0 



8.0 




4.0 6.0 8.0 

Magnetic Field (T) 



10.0 



LLL 



(d) 



A 




J v 



(b) 




1% 



(c) 



Singlet X" 

Triplet X" 

Exciton X 



2meV 



Photoluminescence Energy 



